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SUMMARY 


The  generation  of  thrust  directly  from  electric  power  through  a  gaseous  dis¬ 
charge  is  investigated  with  a  one  dimensional  model.  The  thrust  consists  of  two  parts, 
the  electric  pressure  and  the  electric  wind.  In  terms  of  the  ionization  constants  of 
gases  and  of  the  mobility  constants  of  ions  and  electrons,  both  the  electric  pressure 
and  the  electric  wind  are  obtained  in  simple  dimensionless  form.  At  low  gas  pressures, 
the  thrust  contribution  from  electric  pressure  is  usually  much  larger  than  that  from 
electric  wind.  The  latter  may  be  of  comparable  magnitudes  or  even  much  larger  at  higher 
(atmospheric)  discharge  pressures  and  larger  inter-electrode  distances.  The  overall 
thrust  level  is  proportional  to  the  gas  pressure  in  the  discharge  squared,  p*.  For 
P  -  1  mm  Hg,  the  thrust  is  <  1  dyne/cm*.  The  specific  thrust  per  unit  power  is  no 
greater  than  10  gm  wt/kw.  Thrust  is  derived  from  electric  wind  at  a  much  larger  spec¬ 
ific  power  consumption.  The  low  specific  thruct  alone  rules  out  the  possibility  of 
using  such  a  propulsion  system  as  a  booster.  For  application  in  outer  space,  it  suffers 
from  the  limited  specific  impulse  because  propellant  gas  must  be  carried  aboard. 

There  may,  however,  be  other  interesting  applications,  for  example,  as  a  posi¬ 
tive  control  device  for  hypersonic  vehicles  at  an  altitude  of  100-200  x  103  ft.  A  de¬ 
tailed  feasibility  study  of  this  possibility  has  not  yet  been  made. 
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a  coefficient  in  the  experimental  drift  velocity  expresion;  see  Eqs. (18) (19) » 

1  p.10 

A  saturation  ionization  constant,  cm-1  mm  Hg”1}  see  Eq.(3)»  p.4 

B  ionizing  potential  esu  field/(mm  Hg)s  see  Eq.(3)>  p.4 

d  distance  between  electrodes 

E  electric  field  strength  in  esu,  positive  in  the  direction  of  anode  to  cathode, 

i . e. .  the  positive  x-direction 

j  current  density  of  each  species 

discharge  current  density  esu  current/cm2 

ft 

J  a  dimensionless  discharge  current  density  defined  by  Eq.(20),  p.ll 

p  pressure  of  gas  (plasma)  in  mm  Hg 

r  charge  density  esu/cm3 

Th  reactive  thrust  acting  on  the  system  positive  in  the  negative  x-direction 

Tt  contribution  of  electric  pressure  effect  to  the  reactive  thrust  Thf  see 

Eq.(2l),  p.ll 

Ta  contribution  of  electric  wind  effect  to  the  reactive  thrust  Thj  see  Eq.(22), 

p.  12 

u  mean  velocity  of  the  neutral  plasma  cm/sec  positive  in  the  positive  x-direction 

magnitude  of  drift  velocity  of  particles  relative  to  the  mean  motion  of  plasma 

x  axis  normal  to  the  discharge  electrodes,  anode  at  x  =  0  and  cathode  at 

x  =  d  >  0 

a,p,y,6  ionization  coefficients  defined  in  Eqs. (l) (2) (3) (4) 

y  efficiency  factor  of  momentum  transfer  from  electron  to  neutral  gas»  see 

Eq.(6),  p.5  and  the  following  explanation 

K  same  mean  value  of  v  over  the  entire  discharge  region  defined  by  Eq,(24)»  p.12 

X  external  source  current  density 

w  coefficient  of  cathode  emission  defined  in  Eq.(4)>  p.4 

p  density  of  gas  (in  appropriate  unit) 

£(z)  a  universal  function  defined  by  Eq.(20a),  p.ll 
Subscripts 

-  indicates  quantities  pertaining  to  electrons 

+  indicates  quantities  pertaining  to  ions 

m  indicates  quantities  pertaining  to  metastables 

p  indicates  quantities  pertaining  to  photo-effects 

o  indicates  quantities  evaluated  at  anode  x  =  0 

d  indicates  quantities  evaluated  at  cathode  x  =  d 
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I  PRODUCTION 

When  a  dc  voltage  of  the  order  of  tens  of  kilovolts  is  applied  across  a  pair  of 
parallel  planar  electrodes  constructed  of  light  wire  meshes,  the  electrodes  lift  into 
the  air.1  There  is  little  evidence  of  thermal  expansion  of  the  air  between  the  elec¬ 
trodes.  An  "electro-kinetic"  machine  of  similar  construction2  operating  at  >100  kv  is 
reported  to  turn  a  wheel  with  a  sustained  peripheral  speed  of  the  order  of  10  ft/sec. 
The  novel  feature  of  these  devices  is  that  thrust  is  derived  directly  from  electric 
power  at  room  temperature  without  any  moving  parts  and  that  the  ambient  atmosphere  is 
employed  as  the  propellant.  Might  this  serve  as  an  advanced  propulsion  scheme? 

In  an  electric  field  a  neutral  gas  with  a  slight  degree  of  ionization  (i ,e. ,  a 
neutral  plasma)  is  known  to  move  in  the  direction  of  the  field.  This  is  the  "electric 
wind".5  The  structure  carrying  the  electrodes  will  experience  the  reaction  of  the  gas 
motion,  that  is,  the  thrust.  The  electric  wind  has  been  qualitatively  explained  as 
follows:  A  massive  ion  is  much  more  efficient  than  an  electron  in  transferring  to  the 
colliding  neutral  molecule  the  momentum  and  the  energy  which  the  ion  (or  electron) 
derived  from  the  electric  field  during  the  free  flight  prior  to  collision.  Accordingly, 
the  neutral  molecules  will  gain  a  net  momentum  in  the  direction  of  motion  of  the  ions 
although  the  plasma  may  be  locally  neutral.  The  maximum  effect  of  the  electric  wind 
can  be  estimated  by  ignoring  the  presence  of  electrons. 

However,  if  an  electron  collides  with  a  great  many  neutral  molecules  while 
drifting  back  and  forth  through  a  distance  of  the  order  of  its  free  path,  practically 
all  of  its  momentum  derived  from  the  electric  field  during  the  drifting  motion  may  be 
transferred  to  the  neutral  molecules.  Thus,  when  the  particle  density  is  sufficiently 
large,  depending  upon  the  electric  field  strength,  the  rate  of  momentum  transfer  to  the 
neutral  molecules  per  unit  volume  will  be  E(a+  -  a  ) ,  where  E  is  the  local  electric 
field  strength  and  cr+,c  are  the  magnitudes  of  the  charge  densities  of  ions  and  elec¬ 
trons,  respectively.  Thus,  any  net  force  experienced  by  the  neutral  gas  must  be  due  to 
the  local  excess  of  positive  or  negative  electric  charges.  This  is  the  point  of  view 
of  the  conventional  continuum  analysis  of  plasma  motion. 

Since  the  ionizing  collisions  and  the  drifting  motion  of  the  charged  particles 
are  the  fundamental  elements  of  a  gaseous  discharge,  both  mechanisms  of  thrust  gener¬ 
ation  must  be  present.  Before  analyzing  the  dynamics  of  a  slightly  ionized  gas  under 
electric  discharge  conditions,  a  brief  outline  of  the  physical  processes  and  their 
phenomenological  representations  appears  necessary.  For  a  comprehensive  discussion, 
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the  reader  is  referred  to  standard  textbooks,  in  particular  the  two  volumes  of  the 
Handbuch  der  physik.* 


PHYSICAL  PROCESSES  IN  GAS  DISCHARGE 


Gas  discharge  refers  to  the  passage  of  electric  current  through  a  gaseous  con¬ 
ductor  by  means  of  mass  motion  of  electrons  and  ions  (supposed  to  be  singly  charged 
positively  for  the  sake  of  simplicity)  toward  the  anode  and  the  cathode,  respectively. 
Prior  to  the  breakdown  of  the  gas,  a  steady  discharge  will  have  to  be  maintained  by 
external  means,  for  example,  by  shining  ultraviolet  light  onto  the  cathode  to  cause 
electron  emission.  These  electrons  will  be  accelerated  by  the  electric  field  and  mul¬ 
tiply  through  ionizing  collisions  with  neutral  gas  molecules.  The  current  strength 
will  depend  upon  the  intensity  of  illumination  and  is  usually  very  small,  of  the  order 
of  <  10~1a  amps/cm2. 

When  the  voltage  across  the  electrodes  is  sufficiently  large,  the  discharge 
becomes  self-sustaining*  i . e . ,  a  steady  discharge  current  will  be  maintained  even  if 
the  external  source  of  ionization  is  removed.  From  here  on,  a  slight  increase  in  vol¬ 
tage  will  result  in  a  significant  increase  in  cunent  density.  This  is  the  Townsend 
dark  discharge  ,  where  a  current  density  of  ~10"6  amp/cm2  may  be  reached.  Further  in¬ 
crease  of  the  discharge  current  density  brings  about  the  glow  discharge,  characterized 
by  the  emission  of  visible  l;ght  of  different  patterns.  Attempts  to  increase  the  dis¬ 
charge  current  further  will  bring  about  an  arc  discharge  at  relatively  low  voltages 
and  high  currents  with  hot  cathodes. 

When  the  electric  field  near  the  cathode  is  small,  the  electrons  emitted  from 
the  cathode  collide  with  neutral  molecules  elastically.  When  the  cathode  field  is 
large,  the  electrons  may  gain  a  sufficient  amount  of  energy  during  their  free  flight. 
Their  collision  with  neutral  molecules  may  then  be  inelastic,  resulting  in  many  differ¬ 
ent  products.  The  colliding  molecule  may  be  ionized,  Droducing  an  additional  electron 
and  a  positive  ion.  In  this  manner,  the  number  of  electrons  in  the  gas  increases  with 
distance  from  the  cathode.  The  logarithmic  increment  of  the  electron  concentration 
along  the  distance  from  the  cathode  due  to  the  direct  ionization  by  electrons  is  known 
as  the  first  Townsend  coefficient  a  . 

The  inelastic  collision  of  an  electron  and  a  neutral  molecule  may  lead  to  re¬ 
sults  different  from  ionization.  The  neutral  molecule  (or  the  atom  if  it  is  a  noble 
gas)  may  be  left  electronically  excited  and  may  later  emit  light  by  gcing  to  a  lower 
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state  or  to  the  ground  state.  This  is  the  origin  of  the  visible  glow.  Some  molecules 
(or  atoms)  in  excited  states  may  possess  a  relatively  long  life,  i .  e. ,  metastable 
states,  which  later  may  collide  with  another  neutral  molecule  resulting  in  ionization. 
The  light  emitted  from  the  excited  atoms  may  also  cause  ionization  of  neutral  mole¬ 
cules  through  the  process  of  photo-ionization.  The  rates  of  these  processes  depend 
approximately  linearly  on  the  electron  concentration.  Thus,  the  total  logarithmic  in¬ 
crement  of  the  electron  concentration  a  may  be  taken  as  a  =  a  +  a  +  a  where  a  a i id 

-  m  p  m 

ap  are  the  contributions  from  metastable  atoms  and  photo-ionization,  respectively. 

The  positive  ions,  produced  by  the  previous  ionization  processes,  may  gain  a 
sufficient  amount  of  energy  from  the  electric  field  to  ionize  neutral  particles  upon 
collisions.  The  rate  of  this  ionization  process  will  depend  linearly  on  the  concentra¬ 
tion  of  the  positive  ions,  or  the  ion  current  with  logarithmic  increment  p.  This  pro¬ 
cess  is,  however,  very  inefficient  when  compared  to  ionization  by  electrons.  There  are 
many  other  possible  processes  of  ionization,  including  stepwise  processes,  that  depend 
on  the  degree  of  ionization  to  a  higher  power.  These  higher  order  processes  will  be 
neglected  for  treating  discharge  problems  with  a  small  degree  of  ionization.  Thermal 
ionization,  i . e. ,  ionization  due  to  collision  between  neutral  particles,  and  photo- 
ionization  by  externally  maintained  illumination  are  essentially  independent  of  the 
local  ion  density.  These  may  be  included  in  the  ion  (or  electron)  material  balance  as 
an  additional  source  term.  Hence,  in  terms  of  the  convective  currents,  j  and  j_ 

d)_  dj 

-  -j —  -  aj  +  pjj  +  X  =  +  -t —  (l) 

dx  -  '  J+  ax 

Ions  and  electrons  drift  toward  the  proper  electrode  where  they  are  neutral¬ 
ized.  In  order  that  the  discharge  may  maintain  itself,  the  cathode  must  emit  electrons 
to  maintain  the  subsequent  processes.  The  electron  emission  from  the  cathode  may  be 
produced  by  the  incidence  of  the  positive  ions  or  the  excited  molecules  or  atoms  (meta¬ 
stable  states)  or  of  photons  from  within  and/or  without  the  discharge.  The  electron 
current  density  at  the  cathode  will  hence  depend  linearly  on  the  positive  ion  current 
density  at  the  cathode  with  a  proportionality  constant  y.  The  effects  of  metastable 
atoms  and  photons  are  often  combined  in  the  steady  state  into  one  term  proportional  to 
the  integral  of  the  electron  current  density  over  the  discharge,  (i . ^ ,  proportional 
to  the  rate  of  production  of  photons  and  metastables)  with  proportionality  constant  o. 
The  effect  of  photoerr.i  s  sion  from  the  cathode  due  to  external  illumination  will  be  ne¬ 
glected.  The  thermionic  emission  from  the  cathode  depends  on  the  cathode  temperature, 
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and,  therefore,  the  current  intensity.  We  may  visualize  this  thermionic  emission  as 
contributing  to  y  for  qualitative  purposes.  Hence,  on  the  cathode  surface  (x  =  d) 

rd 

j_(d)  =  yj+(d)  +J  6j_(x)dx  (2) 

o 

All  these  processes  are  probably  present  in  any  gas  discharge  to  some  extent. 
Whether  a  particular  process  may  be  neglected  or  not  depends  on  the  gaseous  pressure 
and  temperature  and  on  the  nature  of  the  electrodes.  For  the  Townsend  discharge,  the 
gaseous  ionization  is  due  primarily  to  electron  avalanche  (a  )  and  the  cathode  emission 
to  positive  ion  incidence  (y) .  For  the  glow  discharge,  the  photo-process  becomes  im¬ 
portant,  and  for  arc  discharge,  the  thermionic  emission  from  the  hot  cathode  predomin¬ 
ates.  Each  kind  of  discharge  will  be  characterized  by  the  predominance  of  some  parti¬ 
cular  process.  The  transition  from  one  type  of  discharge  to  another  is  not  abrupt, 
although  the  transitional  states  may  be  unstable.  The  use  of  the  phenomenological 
constants  is  convenient,  since  the  results  obtained  with  these  constants  are  applicable 
to  a  wide  range  of  situations  so  long  as  the  appropriate  basic  processes  are  identified 
with  these  phenomenological  constants. 

The  ionization  constant  a  for  electrons  are  known  for  many  gases.  a_  must 
be  proportional  to  the  frequency  cf  collision  (~  p  or  p)  and  must  depend  on  the  energy 
that  the  electron  acquires  during  its  free  flight,  1 , c. ,  ~  the  electric  field  times  the 
free  path,  ~  E/p.  Accordingly,  a  semi-empirical  relation  is  often  adopted. 


a_  =  Ap  •  oxp(-  —)  =  a  (3) 

where  A  and  B  are  characteristic  constants  of  a  given  gas.  Such  an  empirical  corre¬ 
lation  gives  a  good  representation  of  the  experimental  data  over  a  fairly  wide  range 
of  p  and  E/p.  The  other  constants  and  are  less  well  known.  Hence  Eq.  (3)  will 
be  adopted  for  a  instead  of  a  with  constants  A  and  B  slightly  adjusted  to  take  and 

a  into  account. 

P 

For  the  processes  on  the  electrodes,  it  is  convenient  to  introduce  a  secondary 
coefficient  w  for  the  cathode  emission  with 


CJ 

u 


y  +  p/o.  +  b/a 


1  +  v  + 


b/  a 


1  -  P/a 


(4) 
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a  quantity  that  can  be  determined  experimentally.  The  magnitude  of  u/a  does  not  vary 
so  critically  with  E/p  as  a  does.  The  major  constituents  of  u/a  for  even  the  simplest 
discharge  processes  are  not  well  understood.  A,  B  and  a/w  are  the  key  parameters  char¬ 
acterizing  the  discharge  process. 

GASDYUAMICS  OF  THE  DISCHARGE 


The  gas  plasma  will  be  assumed  slightly  ionized  and  not  too  much  rarified,  so 
that  the  continuum  equation  of  motion  will  apply.  To  illustrate  the  operating  mechan¬ 
ism,  we  shall  consider  only  the  steady  state  discharge  in  one  spatial  dimension  x  taken 
to  be  normal  to  a  pair  of  parallel  planar  electrodes  with  a  surface  area  large  compared 
with  the  inter-electrode  distance  d.  x  =  0  is  taken  as  the  anode  surface  and 
x  =  d  >  0  the  cathode  surface.  The  negatively  charged  particles  will  be  taken  to  be 
entirely  electrons. 

Under  the  assumption  of  slight  ionization,  the  rate  of  change  of  the  number 
of  neutral  molecules  due  to  ionization,  recombination,  and  diffusion  under  cataphor- 
esis  will  be  neglected  in  the  mass  balance  of  the  neutral  gas.  Hence 


1  -  o 
dx  u 


or  pu  =  const 


(5) 


where  p  is  the  density  of  the  neutral  gas  (or  of  the  plasma  as  a  whole)  and  u  is  the 
drift  velocity  of  the  neutral  gas  with  respect  to  the  electrodes,  positive  in  the  posi¬ 
tive  x  direction  (l . e. ,  from  anode  to  cathode). 

For  the  momentum  balance  of  the  neutral  gas,  we  write 


du  ,  dp 

pu  +  d7  "  E0+  ’  KW- 

=  +  E(a+  -  aj  +  (1  -  *)Ea 


(6) 


On  the  left  hand  side,  the  term  pu  —  stands  for  the  net  out  flux  of  momentum  and  the 
dp 

term  —  the  net  pressure  force.  The  terms  on  the  right  hand  side  stand  for  the  net 
rate  of  momentum  transfer  from  the  charged  particles  to  the  neutral  molecules.  Here 
E  is  the  electric  field  intensity,  positive  throughout  the  region  of  interest  (possible 
exception  in  the  region  of  negative  glow  is  ignored).  C  is  the  magnitude  of  the  charge 
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density  with  subscript  +  and  -  designating  quantities  pertaining  to  ions  and  electrons, 
respectively.  The  momentum  derived  from  the  electric  field  by  the  ions  is  assumed  to 
be  transferred  in  total  to  the  neutral  gas  molecules.  The  momentum  derived  by  the  elec¬ 
trons  is  transferred  to  the  neutral  gas  molecules  with  an  efficiency  vr.  This  efficiency 
ic  will  depend  very  much  on  the  collision  processes,  elastic  or  inelastic,  on  the  par¬ 
ticle  density  of  the  plasma  (pressure),  and  on  the  electron  temperature  which  may  vary 
with  the  distance  from  electrodes.  If  k  =  1,  Eq,(6)  reduces  to  the  conventional  con¬ 
tinuum  momentum  equation  in  which  the  electric  field  exerts  a  force  on  the  neutral  gas 
only  through  the  local  net  charge  density  in  the  plasma.  If  vr  =0,  i  ,e. ,  the  electrons 
are  considered  as  incapable  of  transferring  momentum  to  the  neutral  gas  molecules,  the 
electric  field  will  exert  a  force  on  the  neutral  gas  through  all  the  positive  ions,  but 
not  through  electrons,  as  was  conceived  in  the  mechanism  of  "electric  wind". 

The  convective  currents  j+  and  j  for  ions  and  electrons  in  the  positive  and 
the  negative  x  directions,  respectively,  are 


=*,(vd,-  +  >  0 

j.  -  °-<vd-  -  >  0 


(7) 


where  Vd+  and  V'd_  are  the  magnitudes  of  the  drift  velocities  of  ions  and  electrons, 
respectively.  Vrf+  is  the  mean  drift  velocity  of  ions  toward  the  cathode  relative  to 
the  mean  gas  motion.  Vd  is  that  of  electrons  toward  the  anode.  Both  Vd+  and  Vd_  will 
be  taken  as  known  functions  of  E/p.  The  mean  gas  velocity  u  is  taken  to  be  positive 
in  the  positive  x  direction,  i ,e. ,  toward  the  cathode.  The  field  strength  E  is  always 
positive  in  the  present  sign  convention.  Thus,  the  space  charge  condition  gives 


J_  .  dE  _  j+ _ 

Ljt  dx  +  -  V  -  u  V  ,  +  u 

UT  u- 


(8) 


All  qualities  are  in  electrostatic  units.  The  field  strength,  E  >  0,  decreases  as  x 
increases  from  zero  (proceeding  from  anode  toward  cathode)  because  of  the  abundance  of 
electrons  near  the  anode.  The  cathode  region  is  abundant  in  positive  ions,  where  the 
field  strength  increases  very  rapidly  with  increasing  x  when  the  cathode  is  approached. 

The  conservation  of  ion  and  electron  densities  will  be  written  with  diffusion 
neglected.  It  will  be  clear  later  that  the  neglect  of  diffusion  is  inconsequential  for 
the  present  purpose. 
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^+j+>  =  “j. +  ?h  +  x 

=aj.  +  pj+  +  X 


(9) 


where  ions  and  electrons  are  assumed  to  be  produced  in  pairs  with  ions  all  singly 
charged.  The  meaning  of  a,  £  and  X  are  explained  in  the  previous  sections. 

It  follows  from  (9)  that 

d?(J+  +  O  =0  or  j+  +  j-  =  J0  (10) 

where  Jq  is  the  total  current  density,  independent  of  x.  The  boundary  conditions  at 
the  anode  (x  =  0)  and  cathode  (x  =  d)  are 

0+(O)  =  0  or  j+(0)  =  0 

and 

rd 

j_(d)  =y^+(d)  +  hj_(x)dx  (2) 

J  o 

The  total  voltage  drop  from  the  anode  to  the  cathode,  AV,  will  be 

rd 

AV  =  J  -  E  dx  <  0  (11) 

J  o 

The  electric  quantities  can  be  recognized  as  governed  by  the  three  first  order  differ¬ 
ential  equations  (8)  and  (9)  for  j+(x),  j  (x)  and  E(x)  subject  to  the  three  boundary 
conditions  (2)  and  (ll).  If  we  assume  further 

luA^I,  lu/vd_|  «  1 

and  (12) 

Ap/p  «  1 

to  be  verified  a  posteriori ,  we  note  that  Eqs.(3)  and  (9)  and  boundary  conditions  (2) 
and  (ll)  are  completely  independent  of  the  flow  variables  u  and  p.  Thus,  0+,0  and  E 
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may  be  solved,  as  a  first  approximation  under  conditions  (2)  and  ( 1 1 ) ,  regardless  of 
the  flow  conditions.  The  momentum  balance  of  the  neutral  gas,  Eq.(6),  can  then  be 
solved  when  the  boundary  conditions  of  the  flow  are  specified.  Higher  order  corrections 
to  j+,j  and  E  can  be  evaluated  successively  by  solving  again  (8), (9)  and  so  on.  These 
are  not  necessary  for  our  purposes. 

Combining  (5), (6)  and  (8),  we  can  write 


£(p  +  pu2)  =  +  ^(§£)  +  (1  -  v)  Ep_ 

=  -  *>  E^./(Vd-  -  u) 


(13) 


The  E2/8rr  term  is  the  energy  density  of  the  electric  field  per  unit  volume  which  may 
be  interpreted  as  the  mechanical  force  per  unit  area,  i . e, ,  the  electric  pressure.  The 
next  term  on  the  right  hand  side  of  (13)  represents  the  deficiency  of  momentum  transfer 
by  electrons,  1  ,e. ,  the  electric  wind  effect. 

The  electric  pressure  term  arises  from  the  presence  of  net  electric  charge 
between  the  two  stations.  In  the  vicinity  of  the  anode,  the  electric  field  decreases 
in  the  positive  x-direction  due  to  the  presence  of  excess  electrons.  These  excess 
electrons,  moving  toward  the  anode  (the  negative  x-direction),  transfer  their  momentum 
to  the  neutral  gas  in  the  direction  against  the  flow  of  the  neutral  gas.  The  decrease 
of  electric  field  strength  will,  therefore,  be  accompanied  by  a  decrease  of  the  stream 
thrust  of  the  neutral  gas,  as  indicated  in  Ea.(l3).  Across  the  entire  discharge  re¬ 
gion,  there  will  be  a  net  excess  of  positive  ions,  moving  in  the  same  direction  as  the 
flow  of  the  neutral  gas.  A  net  increase  of  the  electric  field,  accompanied  by  an  in¬ 
crease  of  the  stream  thrust  of  the  neutral  gas,  will,  hence,  result.  The  electric  wind 
term,  being  the  defect  of  the  negative  momentum  transfer  by  electrons, also  acts  in  the 
downstream  direction  to  increase  the  stream  thrust  of  the  neutral  gas. 

Our  object  is  to  evaluate  the  relative  importance  of  the  two  contributions  to 
the  stream  thrust  of  the  neutral  gas.  An  increase  of  the  stream  thrust  of  the  neutral 
gas  means  the  presence  of  a  reactive  thrust  on  the  system  carrying  the  electrodes  per 
unit  sectional  area  of  the  discharge.  This  reactive  thrust  acts  in  the  negative 
x-direction,  i . e, .  in  the  direction  opposite  to  the  flow  of  the  gas.  By  integrating 
Eq . (13)  over  x  =  0  to  x  =  d,  we  obtain  the  following  expression  for  the  magnitude  of 
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the  reactive  thrust 


Th 


Ej_/(Vd_  -  u)]  dx 


(14) 


The  field  strength  E^  at  the  cathode  is  always  higher  than  the  field  Eq  at  the  anode. 
The  efficiency  if  is  never  greater  than  unity.  Hence,  both  the  electric  pressure  and 
the  electric  wind  terms  are  positive,  contributing  to  the  reactive  thrust  on  the  system 
of  electrodes  in  the  direction  from  the  cathode  to  the  anode.  The  net  force  on  the 
plasma  exerted  by  the  system  is  from  the  anode  toward  the  cathode.  In  Fig.  1,  an 
approximate  sketch  of  different  quantities  are  given  without  any  attempt  of  displaying 
the  complicated  fine  details  of  the  discharge  properties. 


ELECTRIC  PRESSURE  VERSUS  ELECTRIC  i/VIND 


The  primary  ionization  constant  a  will  be  taken  as  given  in  Eq.(3)  with  A  and 
B,  given  by  the  experimental  values  in  uniform  fields,  and  E(x)  by  the  local  field 
strength,  p/«,  6/a,  etc.,  will  be  taken  as  some  mean  values  which  are  usually  small. 
Eqs. (9)  can  then  be  readily  integrated 


=(Jo+a)[l 


j  (x)  =  J  e 
-  o 


f(x) 


-  ef(x)]/(l  -  p/a) 

JoP/a+X/a  f(x) 

- - -  [1  -  e 

1  -  P/a 


] 


(15) 


where 

r*  Bp 

f(x)  =  -  (l  -  p/a)  J  a(£))d£  ,•  a  -  Ap  exp(-  -|-)  (16) 

o 

The  condition  at  the  cathode,  x  =  d,  given  as  Eq.(2),  requires  that 


r  <  6/a 

J  [(1  +  Y  +  - T 

0  1  -  p/a 


)  -  (r  ♦  S  ♦  !,] 


(l  +  Y  +  6/a)] 
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Either  when  thermal  ionization  or  photo-ionization  by  external  sources  are  neglected 
(\/a  =  0)  or  when  X/(ay  +  p  +  6)  is  not  »  1,  the  coefficient  of  JQ  must  be  approxi¬ 
mately  zero  since  JQ  in  esu  is  very  large.  Hence 


When  thermionic  emission  is  absent,  the  experimental  values  of  w/o  are  0(10”2  -  10“3) 
varying  slightly  with  E/p.  For  highly  abnormal  glow  discharges  with  large  y,  m/a  may 
be  of  the  order  of  10_1. 

When  a  is  taken  from  Eq,(3)  as  a  known  function  of  E,  Eq.(l5)  and  (16)  express 
j+  and  .1  in  terms  of  E.  When  |u/V^+|  ,  ju/V^  I  «  1  (or  if  u  is  known  from  previous 
iteration)  and  when  V^+  and  are  taken  from  experimentally  known  functions  ol  L/p, 
the  space  charge  equation  (8)  stands  as  a  first  order  differential  equation  for  F(x). 
Thus  E(x)  can  be  solved  in  principle.  The  condition  (ll)  determines  E(x)  completely 
when  the  total  voltage  drop  across  the  electrodes  is  specified.  Solution  of  equations 
in  this  manner,  although  straightforward,  is  tedious  and  does  not  appear  necessary  for 
our  purpose. 

The  drift  velocity  for  ions  way  be  written  according  to  experimental  data5 
(Figure  59-6l)  under  ordinary  vacuum  conditions 


V  =  b+  •  10* (E/p)  for  (E/p)  «  1 

(18) 

=  a+  10*(E/p)*  for  (E/p)  >0(1) 


where  is  expressed  in  cm/sec  and  F./p  in  esu  field/mm  Hg  pressure,  the  coefficients 
b+,a+  are  of  order  of  unity  for  most  gases.  Similarly,  the  drift  velocity  for  elec¬ 
trons  in  gases  can  be  expressed  as 


Vd_  =  b_  109(E/p)^  for  (E/p)  «  1CT1 
=  a  •  10fl(E/p)  for  (E/p)  >  O(lO)-1 


(19) 
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Note  the  difference  in  the  qualitative  dependence  of  and  Vd_  on  E/p. 

While  the  a's  and  b's  are  of  the  same  order  of  magnitude,  the  ratio  Vd+/Vd_  becomes  in¬ 
creasingly  smaller  than  10"  3  with  larger  E/p  (in  esu/mm  Hg).  Thus*  multiply  the  space 
charge  equation  (9)  by  (Vd+  +  u)  •  a  =  Ap  exp[-Bp/E]  *  a+  •  105(E/p)^  and  inte¬ 
grate  from  x  =  0  to  x  =  d.  With  V,./VJ  «  1  and  X/a  «  1  «  J  (in  esu  current),  we 

d+  a-  o 


E, 

r(— ) 

^Bp' 


=J  =  r 


4n- 

A83/2 


Cln(g)  -  lj  rl  -  10”5 


(20)' 


where 


T 


C(2)  =  -  J 

Joo 


1 

T-5/2  e^dr 


(20a) 


£(z)  is  a  universal  function  which  may  be  expressed  in  terms  of  error  function. 
If  V^+  ~  (E/P)s  should  be  adopted  with  s  /  1/2,  Eq.(20)  would  be  modified  as  follows: 
B3/2  -*•  B1+s  and  £(z)  would  be  the  integral  of  T-(2+s^e”T  which  can  be  expressed  in 
terms  of  an  incomplete  P  function  of  order  -(2ts).  Eg. (20a)  therefore  definec  a  univer¬ 
sal  function  £(2)  once  a(E/p)  and  Vrf+(E/p)  of  a  discharge  system  are  selected.  The 

function  £(z)  is  illustrated  in  Fig.  2  for  s  =  1/2.  The  right  hand  side  of  Eq.(20) 

* 

is  a  dimensionless  parameter  J  for  the  reduced  current  density,  JQ/p  ,  involving  the 

characteristic  constants,  a  ,  A,  B  and  a/u,  of  the  gas  in  the  discharge  system. 

+  c 

When  the  dimensionless  field  at  the  anode  -%  and  the  dimensionless  current 
»  P° 

density  J  are  fixed,  Eq.(20)  will  give  the  dimensionless  field  at  the  cathode, 

E^/pB.  Hence,  the  dimensionless  electric  pressure  term  in  (14) 


E  ,  2 


ar^) 

3jt  L  pB' 


(21) 


can  be  calculated. 


This  relation  is,  in  many  respects,  similar  to,  but  different  from,  the  universal  re¬ 
lation  of  Von  Engle  and  Steenbeck6  for  the  cathode  region;  i  . e, ,  V^/p  versus  Jo/pL. 


11 


TP- 33 


The  electric  wind  term  in  (14)  may  be  written  as 


=  J  (1  -  k) 


J  (1  -  k) j_(x)dx 


Let  us  write  j+(x)  =  J  -  j_(x)  in  the  space  charge  equation  (8),  multiply 

through  by  Vd+  +  u  =  a+105(E/p)2  and  integrate  from  x  =  0  to  x  =  d.  We  obtain 

f  (l  -  v)j  (x)dx  ir.  terms  of  J  /p2, (E  ,/p}3/2  and  (E  /p)3/2.  Eq.(20)  may  be  used  to 

J  *  0  0  0 

0  Ed  Eo 

put  J  /p2  in  terms  of  £(up-)  -  C(d  )•  Thus,  the  dimensionless  electric  wind  term  in 

0  bp  Dp 

Eq.(l4)  becomes 

.  P  (l  -  K)dx  , 

T2  Jo  o  +  10'3  rd  (1  -  k)  d  /  E  v 3/2 


B2p2  P2  a_  •  10s 


T  10  f  V  *  ~  *  i  a  /  L  \  ' 

a_  Bi/a  \  6ir  dx  ^Bp 


—  ■  io-«  =-t  A— 1012  _ 

P2  a_  a_  4 n  B1' 2  In 


“  _  i  D^Bp^  '  ^Bp^] 


The  second  part  of  T,/B2p2  is  negligible  compared  with  T-,/B2p2,  This  is  because  the 
integral  in  the  last  term  is  at  most  of  the  same  order  as  T1/B2dz  while  the  coefficient 
is  C(lO"3)  «  1.  The  electric  wind  term  T2/B2p2,  that  may  be  of  any  significance  in 
producing  thrust,  is  the  first  term  in  (23),  which  is  linearly  proportional  to  the  re¬ 
duced  discharge  current  density,  J^/p2. 

E 

Fig.  2  shows  that  is  almost  always  less  than  unity.  The  values  of  A/B2 

for  molecular  and  rare  gases  are  around  10  with  A  in  (cm-mm-Hg)"1  and  B  in  esc  field/mm 
rig.  Except  for  arc  discharges,  a/o  is  substantially  larger  than  1  and  In  q  -  1  is  not 
small.  Consequently,  the  electric  wind  effect  T2  will  contribute  to  the  stream  thrust 
comparable  to  the  electric  pressure  effect  T i  only  when 

J 

p  •  d  •  (l  -  k)  =  p  (l  -  *)dx  >  0(102)  cm  -  mm  Hg  (24) 

J 

0 
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When  condition  (24)  is  fulfilled,  we  have 

+  _ L  C  (E°)2-{  |  Jo  ■  p  ’  d(l  1*1 

B2p2  BJp2  B2p2  Sn-L^Bp'  vBp'  J  p*  a"  •  10® 


(25) 


The  electric  wind  effect  becomes  more  significant  when  the  discharge  pressure 
p  is  higher  and  when  the  inter-electrode  distance  d  is  larger  with  the  reduced  current 
density  JQ/p2  kept  constant.  If  k  «  1,  the  electric  wind  effect  becomes  appreciable 
in  discharges  at  1  mm  Hg  when  the  electrodes  are  several  meters  apart.  At  higher  pres¬ 
sures  and  larger  inter-electrode  distances,  the  electric  wind  effect  may  be  orders  of 
magnitude  larger  than  electric  pressure  when  pd  •  (1  -  k)  •  10“2  »  0(l). 


The  magnitude  of  hr,  i  ,e. ,  the  efficiency  of  momentum  transfer  by  electrons 
upon  colliding  with  neutral  molecules,  varies  considerably,  especially  when  inelastic 
collisions  are  of  importance.  It  is  considerably  larger  than  what  is  indicated  by  the 
ratio  of  electron  mass  to  molecular  mass.  In  a  reduced  electric  field,  as  small  as 
5-10  V/cm  mm  Hg,  the  value  of  the  coefficient  of  energy  transfer  is  0(10_1)  for  elec¬ 
trons  in  molecular  gases  like  N2  because  of  the  tightly  (or  closely)  packed  vibrational 
levels.7  The  corresponding  value  of  the  coefficient  of  momentum  transfer  is  expected 
to  be  of  the  same  order  of  magnitude.  Hence,  k  cannot  be  taken  as  zero  throughout  the 
region  of  discharge  for  the  quantitative  evaluation  of  the  electric  wind  effect  if  the 
electric  wind  effect  should  contribute  appreciably  to  the  total  thrust,  as  indicated 
by  Eq.(25)>  either  for  discharges  at  low  pressures  between  widely  separated  electrodes 
or  for  discharges  at  higher  pressures  (even  sub-atmospheric)  for  moderate  inter-elec¬ 
trode  distances  of  several  centimeters  apart. 

E 

Eq.(20)  is  a  universal  relation  between  the  cathode  field  —  and  the  dimen- 
*  pd 

sionless  current  density  J  .  The  anode  field  to/p8  may  be  recognized  as  a  parameter 
characterizing  the  discharge  in  place  of  AV,  the  overall  voltage  drop  across  the  elec¬ 
trodes.  Thus,  the  thrust  due  to  electric  pressure  (Ti/B2p2)  can  be  given  as  a  set  of 

# 

universal  curves  for  different  values  of  J  and  E  /Bp,  (Fig.  3).  Note  that  the  magni- 

0  # 

tude  of  E  /Bp  is  of  secondary  importance  and  that  when  J  >  1 


(26) 


The  cathode  field  (E^/Bp)  is  also  given  as  a  set  of  universal  curves  (Fig.  4). 

For  estimating  the  orders  of  magnitude,  (l  -  ir)  may  be  taken  as  of  the  order  of  unity 
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so  that 


Ta  ~  A  (p  »  d) 
B2p2  ~  p2  a"  •  108 


(27) 


All  quantities  are  expressed  in  esu  and  ram  Hg  pressure.  These  rough  estimates  will 
guide  the  following  discussion. 

The  characteristic  constants  of  the  gas,  like  A,  B,  etc.,  appear  only  in  the 
dimensionless  parameters.  The  results  of  the  analysis  are  applicable  in  general  for 
any  gas.  We  shall,  however,  take  air  as  a  typical  example  to  facilitate  further  dis¬ 
cussion.  For  discharges  in  air: 

a_  =  1.4  cm  sec"1 (esu  field/mm  Hg) 
a+  =  1.7  cm  sec"1 (esu  field/mm  Hg)1/2 
A  =15  cm-1  mm  Hg"1 
El  =  1.22  esu  field  (mm  Hg)-1 

a/u  may  vary  considerably  even  for  a  glow  discharge,  ranging  from  0(10_1)  to  0(l0“3), 
depending  on  the  cathode  material  and  cathode  surface  conditions.  The  abnormal  glow 
discharge,  in  which  thermionic  emission  from  the  cathode  becomes  appreciable,  may  he 
considered  as  one  with  y  large,  i.e, .  1  >  g  ~  OClO'1).  If  the  thermionic  emission  is 
predominant  at  a  hot  cathode,  as  in  arc  discharge,  there  will  be  severe  variations  of 
gas  temperature  and  pressure  which  are  not  accounted  for  in  the  previous  analysis, 

(p  should  be  replaced  by  the  gas  density  p  in  the  dimensionless  variables.)  We  shall 
not,  however,  seriously  consider  the  arc  discharge  in  the  following  because  its  appli¬ 
cation  to  propulsive  devices  ismore  likely  associated  with  its  thermal  effects. 

The  following  values,  obtained  from  Fig.  3,  illustrate  the  orders  of  magni¬ 
tudes  of  the  quantities  at  i  mm  Hg  pressure  when  the  electric  wind  effect  is  neglected, 
i.e.,  Th  =  Ti . 


Tli  rust 
dynes/cm2 

* 

J 

T1/B2p2 

J 

0 

pA/cm2 

Ed 

V/cm 

-AV 

volts 

Th/ power 

dynes/watt  or  gm  wt/kw 

0.38 

1 .25 

0.25 

9C  - 

-  20 

1000 

0(kv) 

1-10 

1.0 

4.0 

0.8 

300  - 

-  60 

1650 

0(kv) 

1-10 

For  a  discharge  between  electrodes  of  a  given  distance  d,  the  total  voltage 
drop  is  essentially  the  sum  of  the  cathode  fall  (~300  v  for  normal  glow  discharge  in 
air)  and  the  voltage  drop  in  the  positive  column.  The  anode  fall  is  small.  The  reduced 
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Cathode  dark  space  pdn  for  air  is  ~0.2  -0.5  cm-mm-Hg  for  most  electrode  materials. 
Hence*  the  length  of  the  positive  column  is  essentially  the  same  as  the  inter-electrode 
distance  d,  l.e. .  _n  «  1*  The  field  in  the  positive  column  is  smaller  than  the  anode 
field,  but  not  much? different.  A  reasonable  estimate  of  the  voltage  drop  in  the  posi¬ 
tive  column  may  be  Hq  •  d.  Thus,  we  obtain  a  rough  estimate  of  the  overall  voltage 
drop  and  hence  the  power  consumption.  The  calculated  results  agree  with  the  experi¬ 
mental  data  of  Ref.  8,  both  in  the  orders  of  magnitudes  of  numerical  quantities  and  in 
the  qualitative  trends. 


DISCUSSION 

Let  us  first  consider  the  situation  when  the  electric  wind  effect,  TJS  is 
negligibly  small.  The  proportionality  of  the  thrust  T,  to  the  dimensionless  parameter 

#  ft 

J  in  Eq.(26)i  and  the  definition  of  J  in  Eq.(20),  imply  that.: 

1.  A  smaller  discharge  current,  JQ,  will  be  required  to  produce 
a  given  thrust,  Tj,  if  the  ratio  w/a  is  smaller. 

2.  A  smaller  ion  mobility  constant  a, (~larger  ion  mass,  etc.) 
and  a  smaller  saturation  ionization  constant  A  tend  to  reduce 
the  discharge  current  for  a  given  thrust  T1t 

3.  Both  T,  and  Jq  increase  with  the  square  of  the  pressure  of 
the  gas  in  a  given  discharge  system,  i ,e. ,  ~p2. 

The  first  inference  shows  that  abnormal  glow  discharge  with  larger  ut/a,  due  to  ther¬ 
mionic  emission  at  the  cathode,  tends  to  require  more  discharge  current  to  produce 
unit  thrust.  The  abnormal  cathode  fall  also  increases  with  the  discharge  current. 

Hence,  an  abnormal  glow  discharge  will  imply  higher  specific  power  consumption,  although 
a  larger  thrust  may  be  achieved  per  unit  discharge  area.  The  optimum  operating  condi¬ 
tion  will  be  the  normal  glow  discharge  covering  the  entire  electrode.  Larger  thrust 
can  be  achieved  more  economically  by  increasing  the  electrode  area  than  by  operating 
in  the  abnormal  glow  region. 

The  cathode  conditions  (the  field  E,  and  the  current  density  J  )  of  the  glow 

d  o 

discharge  are  not  affected  by  shortening  thp  inter-electrode  distance  within  certain 
limits.  The  entire  column  of  positive  glow  may  be  eliminated  with  the  overall  voltage 
across  the  elpctroaes  reduced  accordingly.  The  power  consumption  in  the  positive  col¬ 
umn  does  not  produce  useful  thrust  when  the  electric  wind  effect  is  negligibly  small. 
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Essentially  the  same  thrust  can  be  derived  from  much  smaller  electric  power  consumption 
by  shortening  the  inter-electrode  distance  to  the  nearly  constricted  condition.  The 
optimum  specific  thrust,  dynes/watt,  is  then  achieved  when  the  overall  voltage  drop  is 
taken  as  the  normal  cathode  fall,  ~3C0  v  for  air.  Thus,  without  the  electric  wind 
effect,  the  thrust  to  power  ratio  cannot  be  expected  to  improve  by  orders  of  magnitude 
from  10  dynes/watt,  or  10  gm  wt/kw. 

The  second  conclusion  indicates  that  for  the  same  thrust,  the  discharge  cur¬ 
rent  density,  and  hence  the  power  consumption,  may  be  reduced  by  using  gases  with  lar¬ 
ger  ion  mass  and  smaller  saturation  ionization  constant  A.  However,  the  range  of 
values  of  the  mobility  constant  a+  and  the  ionization  constants  A  and  B  for  all  gases 
(molecular  or  rare  gases)  are  quite  limited.  Changes  in  orders  of  magnitude  of  the 
specific  power  consumption  cannot  be  expected  by  using  different  gases.  By  proper 
choice  of  electrode  material,  w/a  may  be  reduced  by  several  orders  of  magnitude.  How¬ 
ever,  the  effect  is  felt  only  through  the  logarithm.  Therefore,  reduction  in  discharge 
current  and  power  consumption  by  an  order  of  magnitude  C3n  neither  be  expected  by  se¬ 
lecting  different  electrode  materials. 

The  third  conclusion,  that  increasing  the  discharge  pressure  increases  both 
the  current  density,  Jq,  and  the  thrust,  Th,  in  proportion  to  p2,  means  that  the  speci¬ 
fic  power  consumption  is  not  affected  by  increasing  p. 

Consider  now  the  effect  of  electric  wind  under  the  condition  corresponding  to 
an  electric  pressure  of  1  dyne/cm2  at  I  mm  Hg  and  a  specific  thrust  of  10  dynes/watt. 

If  the  electric  wind  term  T2/B2p2  is  to  contribute  1  dyne/cn2,  we  would  have  (Eq.25), 
p  •  d  =  103  cm-nvrn  Hg.  Corresponding  to  d  ~  10  m,  the  potential  drop  across  the  posi¬ 
tive  column  in  volts  would  then  be 

E+  E+ 

AV+  =  300  B  •  (— )  •  p  •  d  =  365  •  103  •  (^jj) 

where  (E+/pB)  is  the  dimensionless  field  parameter  in  the  positive  column,  which  is 
somewhat  smaller  than  E Q/pb  at  the  anode,  but  is  definitely  of  the  same  order.  Let  us 
take  a  conservative  value  of  E+/pB  =  0.03  (12  v/cm  in  air).  The  voltage  drop  across 
the  positive  column  will  he  10  kv,  i . e, ,  30  times  the  normal  cathode  fall  of  300  v. 

Thus,  the  additional  1  dyne/cm2  thrust  is  obtained  through  an  electric  wind  at  30  times 
the  electric  power  consumption  than  the  1  dyne/cm2  from  electric  pressure.  It  is  thus 
clear  that  the  electric  wind  effect  cannot  improve  the  specific  power  consumption  of 
the  thrust  device. 
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The  thrust  level  that  may  be  derived  from  discharges  at  ambient  condition 
(p  =  760  mm  Hg)  could  be  of  the  order  of  0.5  kg/cm2  if  the  1  dyne/cm2  at  1  mm  Hg  could 

be  extrapolated  according  to  p2  law.  (The  thrust  achievable  from  electric  wind  effect 

is  neglected.)  A  thrust  of  the  order  of  5000  kg  or  10,000  lbs  might  be  obtained  over 

an  electrode  area  of  one  square  meter.  If  the  electric  power  is  supplied  from  a  sep¬ 

arate  source,  it  is  conceivable  that  under  the  condition  of  subnormal  glow  discharge 
(with  little  visible  glow)  corresponding  to  current  densities  of  the  order  of  micro¬ 
amperes/cm2,  a  structure,  weighing  less  than  5  gm/cm2  electrode  area,  may  lift  itself 
against  the  earth's  gravity.  The  demonstration  referred  to  in  the  introduction  is 
clearly  feasible.  This  demonstration,  however,  does  not  imply  the  feasibility  of  such 
a  device  as  a  self-contained  propulsion  system.  If  the  system  is  to  carry  its  own 
power  source  and  lift  itself  from  the  ground,  the  thrust  per  unit  power  must  be  larger 
than  the  specific  weight  of  the  power  plant.  It  does  not  appear  possible  in  the  near 
future,  even  with  excessive  extrapolation  of  the  development  of  present  technology, 
that  there  will  be  any  power  plant  with  specific  weight  of  the  order  of  10  gm/kw. 

•Then  the  discharge  takes  place  in  an  ambient  atmosphere  to  provide  thrust,  no 
propellant  need  be  carriea  aboard  the  vehicle.  The  specific  impulse  is  then  of  no  con¬ 
cern.  The  gas  discharge  system  may  therefore  find  some  interesting  applications  on 
vehicles  for  operation  in  the  lower  atmosphere  with  stand-by  electric  power  supply. 

For  example,  for  vehicles  operating  at  about  105  ft  altitude  (~3C  km)  where  p  ~  10  mm 
Hg,  a  thrust  of  1  kg  wt  can  be  obtained  over  1  square  meter  area  of  control  surface. 
This  may  serve  as  a  simple  positive  control  of  hypersonic  vehicles.  Another  example 
may  be  the  propulsion  of  high  Mach  number  ramjets.  Instead  of  adding  heat  to  the  hot 
air  in  the  diffuser,  just  to  produce  dissociation  of  the  air  rather  than  thrust,  the 
gas  discharge  scheme  will  produce  thrust  without  adding  much  heat  to  the  air.  For  such 
an  application  in  ramjets,  the  inter-electrode  distance  can  be  made  very  large  (if  lib¬ 
eral  supply  of  electric  power  is  available)  to  take  advantage  of  the  electric  wind 
effect.  There  may  be  other  special  circumstances  where  the  gas  discharge  as  a  thrust 
device  may  find  interesting  applications. 

At  altitudes  of  2  •  105  ft  (~60  km)  and  above,  the  atmospheric  pressure  be¬ 
comes  a  small  fraction  of  a  mm  Hg.  Eq.(26),  even  if  applicable,  would  predict  3  thrust 
much  too  small  to  be  of  practical  significance.  When  the  atmospheric  pressure  reaches 
the  range  of  10-3  mm  Hg  (one  micron)  at  about  lfT'  km,  the  free  paths  of  the  ions  and 
electrons  become  so  large  that  they  move  in  an  electric  fieid  as  electron  or  ion  beams, 
rather  than  drift  as  swarms  of  random  particles.  They  may  also  attain  energy  levels 
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high  enough  to  cause  many  new  phenomena.  Therefore,  discharges  through  the  highly  rar¬ 
efied  "atmosphere"  at  such  altitudes  cannot  be  discussed  within  the  present  context. 

If  the  gas  discharge  in  the  ordinary  sense  is  to  be  applied  at  such  high  altitudes  or 
in  the  environment  of  outer  space  as  a  propulsion  device,  the  gas  must  be  stored 
aboard  the  vehicle  under  a  pressurized  chamber  and  allowed  to  flow  through  the  region 
of  discharge  into  outer  space.  The  electric  pressure  is  likely  to  be  much  smaller  than 
the  gas  storage  pressure.  The  electric  wind  effect  can  neither  be  counted  on  to  in¬ 
crease  greatly  the  thrust  level,  because  the  loss  of  neutral  gas  through  diffusion  into 
outer  space  -vi 1 1  be  excessive.  The  low  specific  impulse  will  put  the  gas  discharge 
system  in  a  rather  unfavorable  competitive  position  for  space  applications. 

In  conclusion,  the  glow  discharge  in  gases  (without  resorting  to  thermal  ef¬ 
fects  like  the  arc  discharge)  is  limited  in  its  potentiality  as  an  advanced  propulsion 
device.  Within  a  reasonable  expectation  cf  the  specific  power-weight  ratio  of  electric 
power  plants,  the  gas  discharge  system  cannot  serve  the  purpose  of  a  primary  propulsion 
system.  However,  there  may  be  other  interesting  applications,  like  the  auxiliary  thrust 
device  for  vehicle  control,  and  the  propulsive  device  for  high  Mach  number  ramjets. 
Detailed  study  for  the  specific  function  has  yet  to  bp  made 
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